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Abstract 

We report a parton distribution function (PDF) analysis of a complete set of hadron scattering 



i-C ' data, in which a color-octet fermion (such as a gluino of supersymmetry) is incorporated as an extra 

■ parton constituent along with the usual standard model constituents. The data set includes the 

most up-to-date results from deep-inelastic scattering and from jet production in hadron collisions. 



Another feature is the inclusion in the fit of data from determinations of the strong coupling a s (Q) 
at large and small values of the hard scale Q. Our motivation is to determine the extent to which 



(f-) , the global PDF analysis may provide constraints on the new fermion, as a function of its mass 

and a s (Mz), independent of assumptions such as the mechanism of gluino decays. Based on this 
analysis, we find that gluino masses as low as 30 to 50 GeV may be compatible with the current 
hadronic data. Gluino masses below 15 GeV (25 GeV) are excluded if a s (Mz) varies freely (is 
equal to 0.118). At the outset, stronger constraints had been anticipated from jet production 
cross sections, but experimental systematic uncertainties, particularly in normalization, reduce the 
discriminating power of these data. 

PACS numbers: 13.87.-a, 14.80.Ly, 13.85.Qk 
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I. INTRODUCTION 



Heavy color-octet particles are postulated in theories of beyond-the-standard-model 



(BSM) phenomena, including supersymmetry (SUSY) universal extra dimensions |2j, 
Randall-Sundrum [3[, and Little Higgs models |4J. Direct searches for such states are usually 
guided by aspects of the production and decay dynamics in the particular BSM approach. 
Analyses of search data have so far produced various bounds on the masses of the states, 
often conditioned by model-dependent assumptions {5 15]. Different constraints, such as the 



SUSY gluino mass bounds m- q > 26.9 GeV 



16J and 51 GeV [17| at 95% confidence level 



(C.L.), are based on the analysis of LEP event shapes in soft-collinear effective theory and 
other quantum chromodynamics (QCD) resummation formalisms. Constraints such as these 
may depend on theoretical modeling of nonperturbative hadronization and the matching of 
hard-scattering and resummed contributions, similar to the determination of a s (Mz) from 



LEP data in QCD 



18 



231 ] . we examine the possibility that a 



3H22J. In a previous publication 
global analysis of hadron data, within the framework of parton distribution function (PDF) 
determinations, can be used to derive constraints on the existence and masses of color-octet 
fermions, independently of other information on such states. Global analysis has discriminat- 
ing power for several reasons: one is that new colored states modify the evolution with hard 
scale Q of the strong coupling strength a s (Q). Second, in perturbative QCD, the coupling of 
a color-octet fermion to quarks and gluons alters the set of evolution equations that governs 
the behavior of all parton distribution functions, thus affecting many hadron scattering cross 
sections. Moreover, production of the color-octet states will affect relevant observables, such 
as jet rates, whose cross sections are included in the global fits. 

The specific case of a gluino from supersymmetry is included as an extra degree of freedom 
in our earlier work j^. We refer to the PDFs obtained in that publication as "SUSY PDFs", 
although our analysis is applicable to a broader class of standard model (SM) extensions. 
In Ref. 231 ] . a lower bound on the gluino mass rrtg is obtained in terms of an assumed 
value of a s (Mz) at Z boson mass Mz- For the then standard model world-average value 
of a s [Mz) — 0.118, gluinos lighter than 12 GeV were shown to be disfavored, whereas the 
lower bound was relaxed to less than 10 GeV (less than 2 GeV) when a 8 {Mz) was increased 
above 0.120 (0.127). 
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In this paper, we use new hadron scattering data incorporated in the next-to-leading order 
(NLO) CT10 general-purpose PDF analysis 24|, along with a new approach for incorporating 
the variation of a s (Q) into PDF determinations 25j, to obtain improved bounds on the mass 
of a relatively light gluino. The essential new elements are these: 



New Tevatron jet data [26M28I] and combined DIS data 29] from HERA. In a 

global QCD analysis, the presence of light gluinos is revealed primarily by modifications 
of a s (Mz), the gluon PDFs, and the charm and bottom quark PDFs, generated ra- 
diatively above the respective heavy-quark thresholds. We include the latest hadronic 
scattering data sensitive to such modifications. The most stringent constraints on 
the gluon PDF are imposed by electron-proton deep inelastic scattering (DIS) data 
at x < 0.1 and single-inclusive jet production data from the Tevatron pp collider at 
x > 0.1. The study reported here incorporates up-to-date information from the com- 
bined HI Collaboration and ZEUS Collaboration data on deep inelastic scattering 
at HERA-1 29j, as well as single-inclusive jet data from the Tevatron Run-II analy- 



ses 126 



28] . Hard scattering contributions of massive gluinos, with full dependence on 



the gluino's mass, are included in the jet production cross sections we use, the only 
process we examine where these contributions are large enough to be relevant at NLO 
accuracy. 

Floating a s (Mz). Our fits are performed by treating a s (Mz) at the mass Mz of 
the Z boson as a variable parameter of the standard model. We constrain a s [Mz) by 
requiring that the fitted a s (Q) agree with its direct determinations at low energy scales 
(Q < 10 GeV) and at Q = Mz, within the quoted uncertainties of these measurements. 
Virtual gluino contributions result in a slower evolution of the QCD coupling strength 
a s (Q) at scales Q above the gluino mass threshold. By including data that constrain 
a s at low and high Q scales, we effectively probe for deviations from pure QCD. We 
find, in particular, that the value of a s [Mz) = 0.123 ± 0.004 derived in some analyses 



of LEP event shapes 



201 ] can be accommodated if gluinos have mass of about 50 GeV. 



The remainder of this paper is organized as follows. In Sec. [TT1 we describe the role of new 
color-octet fermions in a global QCD analysis. The incorporation of data on a s (Q) within 
the global fits is discussed in Sec. 1111} where we also present the values of a s {Q) at high and 
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low Q used in our fits. Our simultaneous global fit to hadronic scattering data and a s (Q) 
is described in Sec. HVt where we also examine the effects of an additional gluino degree of 
freedom on the PDFs. We present figures that show the relative magnitudes of the PDFs 
and the variation of their momentum fractions with gluino mass and hard scale. 

Section |V] contains the results of our detailed comparison with data. We present figures 
that show the variation of the values of x 2 in the global analyses, as a function of gluino mass, 
for both floating and fixed a s {Mz)- Section [V] also includes the comparison of our calculated 
cross sections with jet data from the Tevatron collider and a discussion of the systematic 
uncertainties that limit the constraining power of these data. The sensitivity of jet cross 
sections at the LHC to the presence of gluinos is examined in Sec. IVII Our conclusions are 
presented in Sec. IVII I The appendices contain an analytic expression for the evolution of 
the strong coupling a s (Q) in terms of the SM and SUSY degrees of freedom; expressions for 
the contributions of massive gluinos to the jet production cross sections; and parton-parton 
luminosity functions for various combinations of SM partons and gluinos. 

Based on our analysis, we conclude that gluino masses as low as 30 to 50 GeV may 
be compatible with the current hadronic data, depending on the value of a s [Mz)- For a 
floating a s (Mz), gluinos lighter than 15 GeV are excluded. For an assumed fixed value 
a s (Mz) = 0.118, the world average value used in many phenomenological analyses, gluinos 
lighter than 25 GeV are disfavored. 

We acknowledge that a gluino as light as ~ 50 GeV is not typical in phenomenological 
models of SUSY breaking, nor of the results of experimental direct search analyses based on 
specific models of SUSY breaking and assumptions about mass relationships among SUSY 
states [15j. As long as the SUSY neutralino x° is lighter than the gluino, the typical decay 
process for a light gluino is g — > qqx°, where q stands for a SM quark. Missing energy 
would signal the presence of a neutralino. However, for a small mass splitting mg — m^o, the 
gluino's decay into missing energy and soft quark jets would be undetected. The analysis 
reported here is complementary to other approaches for bounding the gluino mass, and it is 
in some respects more general in that we make no assumptions about the gluino decay. 

Precise determination of a s (Mz) and proton PDFs are essential ingredients for obtaining 
reliable predictions from perturbative QCD calculations. Such calculations are key for the 
general physics program and for new physics searches at the CERN Large Hadron Collider 
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(LHC) and Fermilab Tevatron collider. As we show here, these ingredients themselves may 
be affected by non-SM contributions, at all values of the momentum fraction x, as a result 
of the global interconnections in PDF analyses. The determination of the QCD coupling 
a s and of the gluon PDF from the Tevatron or LHC single-inclusive jet data, such as in 

I30I ]. may be sensitive to scattering of color- 



recent Tevatron Run-2 measurements 



octet fermions in the ways discussed in Sec. IVI1 As a result of our work, we determine new 
sets of PDFs that include a relatively light gluino as a hadron constituent. 



II. COLOR-OCTET FERMIONS IN A GLOBAL QCD ANALYSIS 

Under well-defined conditions, a relatively light strongly-interacting fundamental particle 
may be treated as a constituent of the colliding hadrons. It will share the momentum of 
the parent hadron with the standard model quark, antiquark, and gluon partners. The 
experimental consequences of this picture become evident when the parent hadron is probed 
at a sufficiently large hard scale. For example, the charm quark c and bottom quark b are 
treated appropriately as partonic constituents of hadrons when the characteristic energy 
scale Q exceeds the mass of the heavy quark m q . Likewise, when Q greatly exceeds the mass 
of a new strongly-interacting particle, this object must also be incorporated as a hadronic 



constituent. We refer to Ref. 



231 ] for an exposition of the PDF analysis in which a gluino is 



included as an additional partonic degree of freedom. 



As in Ref. 



23] , we take the gluino as the only colored non-SM degree of freedom that needs 



to be considered. In some models of SUSY breaking, such as split supersymmetry [3l|, |32j, 
the squarks are much heavier than the gluinos, and therefore could be omitted from our PDF 
analysis. Moreover, as illustrated in Eq. (A3) of Appendix A, color-octet spin- 1/2 fermions 
have a greater impact on the evolution of the strong coupling a s than color triplet scalars, 
such as squarks. 1 

The presence of a light gluino g modifies the PDF global analysis in three ways. 

1. The gluino changes the evolution of the strong coupling strength a s (Q), as the scale Q 
is varied. This influence is implemented in our results, and we provide details on the 



1 While bottom squarks (b) can be relatively light in some models 33( , their contribution to DIS and other 
relevant cross sections can be neglected, cf. Ref. Q. 
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running of oe s (Q) in Appendix [A] The constraints on the gluino mass from our global 
analysis depend significantly on the value of the strong coupling strength a s (Mz). 

2. The gluino provides an additional partonic degree of freedom that shares in the nu- 
cleon's momentum. It alters the coupled set of Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi equations that govern the evolution of the parton distributions, 



d 



dQ 1 



y 



g(x,Q) 

\ g(x,Q) J 

f P^°(x/y) P£ g LO (x/y) P^{ x /y)\ ( 



a s (Q) 
2vr 



NLO, 



x/y) P g N g LO (x/y) P g f{x/y) 
{ P${x/y) P g L g °{x/y) PiP {x/y) ) 



S(y,Q) 
g{y,Q) 
\ g{y,Q) J 



i=u,d,s, 



Here H(x,Q), g(x,Q), and g(x,Q) are the singlet quark, gluon, and gluino distribu- 
te, Q) are the quark and antiquark distributions for 



tions, respectively; qi(x,Q) and cji 
a flavor i. The previous analysis 23| shows that the gluino's contribution is small in 
the momentum fraction range x > 10~ 5 , g(x,Q) g(x,Q), and g(x,Q) q(x,Q). 
NLO variations in the relevant SUSY cross sections are small and comparable in size 
to variations associated with next-to-next-to-leading order (NNLO) SM contributions. 
Therefore, the leading-order (LO) approximation for the splitting functions and hard 
scattering amplitudes of SUSY terms is numerically adequate, when combined with 
NLO expressions for SM contributions. 

3. At energies above its mass threshold, a color-octet fermion contributes to hard scatter- 
ing processes as an incident parton and/or as a produced particle. However, as argued 



in Ref. 



23] , in the absence of light squarks, gluino hard-scattering contributions to 



DIS and Drell-Yan process are of next-to-next-to-leading order and negligible in the 
current study. At the same time, the hard-scattering gluino terms contribute at the 
LO in single-inclusive jet production, so that it is essential that we include the gluino 
in the corresponding hard scattering matrix elements of jet cross sections. 
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FIG. 1: LO scattering diagrams with gluinos in the initial or final state. The double lines stand 
for the squark exchange contributions that we neglect in our approximation. 



The 2 — y 2 hard scattering contributions with two gluinos in the initial or final states 
are illustrated in Fig. [TJ We assume that the masses of the squarks are large enough 
that diagrams containing a squark propagator are negligible. The remaining SUSY 
diagrams can be evaluated in the the S-ACOT factorization scheme 



35| . in order 



to simplify treatment of the gluino mass dependence. In this scheme, gluino mass terms 
are retained in diagrams with two final-state gluinos in the subprocesses gg —y gg and 
qq — y gg. Explicit scattering amplitudes in these channels are documented in Eqs. ( 1B2I) 
and ( 1B3I) of Appendix [B] Massless amplitudes are used for the remaining 2 — y 2 hard 
scattering subprocesses, in which one or two gluinos are present in the initial state, 
and whose contributions are proportional to the gluino PDF g(x, Q). This arrangement 
captures the full gluino mass dependence, while including the mass terms only in the 
essential scattering amplitudes. 



III. QCD COUPLING STRENGTH AS A FITTING PARAMETER 

Since the range of m- g values allowed by the global fits depends strongly on the assumed 
value of a s (Mz), we do a simultaneous fit to hadronic data and to data on a s (Q) in this 
work. A judicious choice is required therefore of the set of data on a s (Q). 

Our approach is to fit the global set of data using a s (Q) as a floating parameter, con- 
straining it with additional data on a s (Q) measurements at Q < 10 GeV (i.e., in the range 
where gluino contributions are excluded by the previous analysis), and at Q = Mz (in e + e 



hadroproduction at LEP). This approach is similar to the floating a s (Mz) fit in Ref. 251 ]. 



However, we constrain at s (Q) at two distinct Q values, to probe for deviations of its running 



from the SM prediction. 



A. Low-energy constraints 

The QCD coupling constraint at low Q = 5 GeV, 

a s (Q = 5 GeV) = 0.213 ±0.002, (2) 

is obtained as a weighted average of three precise determinations of a s at comparable energies: 

a s (Q = 5 GeV) = 0.219 ± 0.006 from r decays, (3) 
a s (Q = 5 GeV) = 0.214 ± 0.003 from heavy quarkonia, (4) 
a s (Q = 5 GeV) = 0.209 ± 0.004 from lattice QCD. (5) 

These values are reconstructed by QCD evolution to the common scale Q = 5 GeV of the 
published a s values provided at different energy scales, 

(a s ) T = 0.330 ± 0.014 at m T = 1.77 GeV, (6) 
(a s ) QQ = 0.1923 ± 0.0024 at M qQ = 7.5 GeV, (7) 
Kkttkc = 0.1170 ± 0.0012 at M z = 91.18 GeV. (8) 



The value of (a s ) T is determined from measurements of r decays 36|; (o: s )qq comes from 
heavy quarkonium decays [^j; and (a s )i a ttice is obtained from lattice computations 37]. 

The r decay and heavy- quarkonium determinations of a s can be reasonably assumed to 
be independent of gluino effects. Even if very light gluinos (~ 10 GeV) were present, the 
value of a s in these measurements would not be affected. The lattice QCD value (tto)iattice 



is also determined at Q < 10 GeV from the energy levels of heavy quarkonia [37J], and 
then evolved by the authors to Q = Mz assuming the SM /3-function. We reconstruct the 
"directly measured" lattice QCD value at Q = 5 GeV (independent of the gluino effects) by 
backward SM evolution. We then combine the lattice QCD value with the other two low-Q 
measurements, evolved to the same scale using the SM /3-function, to obtain a composite 
data input to the fit. 
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B. Z pole constraints 



If nig < Mz, the value of a s {Mz) extracted from the LEP e + e~ hadroproduction data 
could differ from the value obtained from SM fits. On the other hand, various determinations 
of a s (Mz) from Z boson width and hadronic event shapes [l8 - 



22 



38 



show no obvious 



need for BSM contributions. Thus, if gluinos are lighter than Z bosons, their contributions 
to the LEP observables are of the order of theoretical uncertainties from other sources. 
Notably related to assumptions about nonperturbative hadronization in LEP observables, 
these uncertainties remain substantial and produce central values of a s [Mz) ranging from 



0.1135 



21 



22[ to 0.1224 [20j. To deal with this issue of choice, one solution is to include 
available values of a s (Mz) derived from the Z width and/or event shape measurements, 
assuming that gluino contributions for these measurements are comparable with the current 
experimental plus theoretical uncertainties. 

Absent a gluino lighter than the Z boson (i.e., if only SM particles contribute at Q < Mz), 
NLO evolution of the composite low-Q value in Eq. to the Z pole results in a s [Mz) 



close to 0.118. Global analysis of hadronic scattering alone a 



so leads to a preferred value 



a s (M z ) = 0.118 ± 0.005 at 90% C.L., cf. recent CTEQ fits [25]. 

If the gluino is lighter than Mz, the resulting evolved value at Q = Mz is higher. For 
example, the evolved a s (Mz) is 0.126 or 0.121, if m g is 20 or 50 GeV. This variation is 
illustrated in Fig. [21 showing the dependence of a s (Q) on the scale Q in the absence of light 
gluinos (solid line) and with gluinos of mass m g = 50, 25, 10 and 5 GeV. In the figure, we 
show the low-Q constraint (the left data point), as well as one of available constraints at the 

n 

Z pole, a s {Mz) = 0.123 ± 0.004 [20J. As seen in the figure, a light gluino with a mass of 
m g = 10 GeV cannot simultaneously accommodate the low-Q and high-Q constraints. On 
the other hand, gluinos with mass about 50 GeV are compatible with both constraints, and 
are even preferred if the high-Q constraint on a s is larger than 0.118. 

To illustrate typical possibilities, we therefore present two kinds of fits in this paper: one 
in which a fixed value of a s [Mz) = 0.118 is assumed; and the other in which a s (Mz) varies 
and is constrained by an assumed high-Q data point, 

a s (M z ) = 0.123 ±0.004, (9) 



compatible with [20 ]. 



a" 




Q [GeV] 

FIG. 2: Running of the strong coupling as a function of the scale Q. The red solid line represents 
the SM evolution, while the dashed lines are plotted for rrig = 50,25,10,5 GeV. The points with 
the error bands represent the low-Q and high-Q constraints, given in Eqs. (T2J) and ([9]), respectively. 

C. Log-likelihood function for coupling strength constraints 

With the additional constraints on the running coupling, the total log-likelihood function 



Xtot is 



2 _ 2 |_ 2 
Xtot Xh.s. * Xa s : 



(10) 



where x\ s * s the X 2 contribution of the hadron scattering (h.s.) experiments, i.e., DIS, 
vector boson production, and jet production; and xi s * s the contribution from the direct 
constraints on a s : 



x 2 

Aa s 



OLs \exp — OLs \th 



5a 



(0| 



11) 

i=l \ uOL s \exp 

In this equation, N as is the number of data points constraining a s ; N ag = 2 in our case. 
&s^\exp an d SaLslexp are the central value and error of the experimental measurements in 
Eqs. (j2J) and ([9]); ai \th are the respective two-loop theoretical values. We assume that 
an increase in \ 2 by 100 units above the best fit value corresponds to approximately 90% 



CL. error, in accordance with the convention of the previous CTEQ6 analysis 



and 2004 



gluino study |23| . To match this convention, the a s contribution \ 2 a ls included with a factor 



10 



A = 37.7, so that a deviation of a® \th by l.66a^\ exp (90% C.L.) corresponds to A X l s » 100. 

IV. GLOBAL FITS 

In this section we describe our simultaneous global fit to hadronic scattering data and 
a s (Q), and we examine the effects of an additional gluino degree of freedom on the PDFs. 
Our SUSY fits include the same set of data as the latest CT10 fit of parton distribu- 



tions 24j. A total of 2753 data points from 35 experiments is included. Besides the data 



studied in the previous CTEQ6.6 analysis 40], the new analysis includes the combined DIS 



data from HERA-1 



2n 



and single-inclusive jet data from the Tevatron Run-2 analyses [26- 
The new data provide important constraints on the gluon PDF, the parton density that 
is most affected by the gluinos. The charm and bottom PDFs are also affected, since they are 
generated by DGLAP evolution from the gluon PDF above the initial scale Qo = m c = 1.3 
GeV. 

The ratios of the best-fit gluon and charm PDFs in the SUSY sets to their counterparts 
in the standard model CT10 set, fsusY(x,Q)/fcTio{x,Q), are shown as dashed curves in 
Figs. [3] and HJ at Q = 2 and 85 GeV, for two values of the gluino mass, m- g = 20 and 50 
GeV. The normalized CT10 uncertainty bands are shown also, defined as 

fcTio(x, Q) ± 5±f C Tio(x, Q) 



fcrio(x,Q) 



(12) 



in terms of asymmetric PDF uncertainties 5±fcTio{x,Q) [2J, |4l|. Figure [3] pertains to fits 
with a floating a s (Mz), whereas Fig. H]is based on a fixed a s (Mz) = 0.118. 

If a s [Mz) varies (Fig. [3]), modifications in the gluon distribution are moderate at most. 
Some differences with the CT10 predictions are observed at large x, notably in the range 
x > 0.01 in g(x, Q) at Q = 85 GeV and in c(x, Q) at Q = 2 GeV. The difference is larger for 
a lighter gluino with rrtg = 20 GeV. Other PDFs exhibit smaller differences, all contained in 
the standard model uncertainty band. 

For a fixed a s (Mz) = 0.118 (Fig. HJ), the differences with CT10 are substantial. At 
Q =2 GeV, the SUSY PDFs lie outside of CT10 error bands for x as low as 10~ 3 . At Q =85 
GeV, the difference persists at x > 0.01 — 0.05. Large differences between the SUSY and 
CT10 PDF's in the case of a fixed strong coupling are attributed to sizable deviations from 
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g(x,Q) Q = 2GeV, CT10: red error band 
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FIG. 3: Ratios of g(x, Q) (upper row) and c(x, Q) (lower row) distributions in SUSY fits with 
floating a s (M z ) and the CT10 fit at Q = 2 GeV (left) and Q = 85 GeV (right), for the gluino mass 



m-g of 20 and 50 GeV. 



SM running of a s and compensating adjustments in g(x, Q) observed for relatively light 



gluinos; cf. Figs. 12b and 5 in Ref. 



23] 



The effect of the gluino on the standard model quark and gluon PDFs can be significant, 
even if m g is large compared to m c and m^. 2 Because the gluino is an active constituent 
of the proton, it carries a finite momentum fraction, taken from the other non-SUSY par- 
tons, primarily the gluon. This feature is evident in Table |T] where we display the partonic 
momentum fractions for gluino masses m g = {20, 50, 100} GeV. 

Gluinos draw most of their momentum fraction from the gluon, since the primary coupling 
is via the process g — > gg. The influence on the quarks is a second-order effect transmitted 



2 We use m c = 1.3 GeV and m& = 4.5 GeV. The up, down and strange quark masses {m u , nid, m s ) do not 
play a role in the evolution, as they are less than the initial evolution scale Qo = 1.3 GeV. 
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g(x,Q) Q = 2 GeV, CT10: red error band. a s (M z )=0. 1 IS 
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FIG. 4: Same as Fig. but for a fixed a s (M z ) = 0.118. 



Momentum fractions for Q = 100 GeV in percent 


m g [GeV] 
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b 


20 


2.8 


3.9 


3.4 


44.3 


21.8 


11.4 


3.0 


1.8 


1.2 


50 


1.2 


3.9 


3.4 


45.8 


21.8 


11.4 


3.1 


1.9 


1.2 


100 





3.9 


3.4 


47.1 


21.7 


11.4 


3.0 


1.9 


1.2 



TABLE I: Momentum fraction = J dxx fi(x, Q) for each partonic flavor i at scale Q = 100 GeV. 
Momentum fractions for {s, c, b} are not shown and must be included to satisfy the sum rule. 

through the gluon. At Q = 100 GeV, the momentum fraction of the lighter gluinos (nig ~ 
20 GeV) is comparable to that of the strange quark, even though the gluino mass is an order 
of magnitude larger. For m g ~ 50 GeV, the momentum fraction of the gluino is comparable 
to that of the bottom quark. The magnified impact of the gluino on the QCD evolution, 
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m~=50GeV, Q=100GeV 




FIG. 5: PDFs for various flavors at Q=100 GeV, for m g = 50 GeV. 



compared to the usual quark flavors, can be understood from a comparison of the g — >■ g 
splitting kernel, 

Pg-+g(x) =3[(1-X) 2 +X 2 ], (13) 

with the usual gluon-quark splitting function 

P g ^ q {x) = - xf + x% (14) 

The effect of the gluino as a hadronic constituent in the QCD evolution is thus equivalent 
to that of 6 quark flavors, P g -> 9 = 6 P g ^ q - 

As an illustration of the relative magnitude of the gluino PDF, Fig. [5] displays PDFs for 
various parton flavors as a function of x for our m g = 50 GeV PDF set and a hard scale of 
Q = 100 GeV. At x > 0.001, the gluino PDF is about equal to the bottom quark PDF, the 
smallest of the quark PDFs. For smaller x, it grows in magnitude and catches up with the 
other quark PDFs at x = 10~ 5 , as a consequence of its faster DGLAP evolution. Parton- 
parton luminosities dependent on the gluino PDF, useful for computations of cross sections, 
are plotted in Appendix O 
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SUSY analysis with a fixed a s (M z ) =0.118 


m s [GeV] 


Xh.s. 


Xtot 


X 2 /npt: HERA-1 


X 2 / n P t: jet prod. 


a s (M z ) 


10 


3154 


12550 


1.31 


1.24 


0.118 


20 


3030 


7882 


1.24 


1.19 


0.118 


50 


2923 


3788 


1.18 


1.10 


0.118 


00 


2918 


3004 


1.16 


1.09 


0.118 


SUSY analysis with a floating a s (Mz) 


m-g [GcV] 


Xh.s. 


Xtot 


X 2 /npt: HERA-1 


X 2 /npt: jet prod. 


a s (M z ) 


10 


2892 


3124 


1.14 


1.06 


0.132 


20 


2897 


2958 


1.15 


1.06 


0.127 


50 


2896 


2901 


1.15 


1.03 


0.121 


oo 


2918 


2960 


1.16 


1.09 


0.118 



TABLE II: \ 2 values in the global analyses with a floating and fixed a s (Mz), for various gluino 
mass values. 



V. COMPARISON OF THEORY AND DATA 



In this section, we show the results of our global fits, the constraints we obtain on the 
mass of a gluino, and the impact of a gluino degree of freedom on the analysis of jet data. 

The figures in the previous section show that the SM+SUSY PDFs disagree with CT10 
PDFs if gluinos are lighter than 20 GeV, indicating that the SM+SUSY PDFs for these 
gluino masses cannot describe the global hadronic data well. Gluinos with somewhat larger 
masses can be accommodated, or may be slightly preferred to the pure SM case, depending 
on the value of a s [Mz)- These points are illustrated in a different way by the summary of 
values of x 2 i n Table HT| for m- g =10, 20, and 50 GeV, as well as for the standard model case 
(equivalent to rrig = oo). The table shows the log-likelihood values Xh.s. an< ^ Xtoti without 
and with the imposition of a s constraints, as defined in Eqs. ( fTOl) and (TTT1) ; as well as x 2 P er 
number of data points for HERA-1 DIS 29( and Tevatron Run-1 and Run-2 single-inclusive 



jet cross sections 



26-28 



42 



431 ] . In the fit with a floating a s , the best-fit a s [Mz) is also 



shown. A comparison of the upper and lower halves of the table shows that the relation 
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SUSY fits with a fixed a s (M z )=0.1 18 




gluino decoupling - 



10 20 30 50 

m ? [GeV] 



100 



200 



FIG. 6: Values of A%? vs. m g are shown for a fixed value of a s (Mz) = 0.118 for the 2004 study 
(blue dashed line) and our new one (red solid line). 



between x 2 an d m g depends on whether a s (Mz) is fixed or floating. 

Fixed a s (Mz). In a fit with a fixed a 8 (Mz), only constraints from the hadronic data, 
associated with the term \\ s ( an d not with the total Xtot) P ia Y a meaningful role. The upper 
half of Table HT1 shows x 2 values from a fit with fixed a s (Mz) — 0.118. 3 In this case, the 
gluino's effect of slowing the evolution of a s (Q) from Q = M z to Q = 5 GeV runs into strong 
disagreement with the low-Q constraint; Xtot g rows quickly as m g decreases, corresponding 
to a difference of many standard deviations between the measured and predicted a s values at 
Q = 5 GeV. More importantly, the hadronic data by themselves disfavor very light gluinos, 
with m g = 25 GeV or less excluded according to the criterion A% 2 = X 2 susY( m g) ~~ Xctio < 
100 applied to x\. s .- 

Floating a s (M z ). The values in the lower half of Table E] are for SM+SUSY fits with 

3 This value, compatible with the current world average, is about la below a s (Mz) = 0.123 ± 0.004; hence, 
the SM fit with this a s {Mz) value has a higher xiot (i n the last line of the upper table) than a fit with a 
floating a s (Mz) (in the last line of the lower table). 
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SUSY fits with a floating 0C s (M z ) 



200 




2010 study 
2004 study 



150 - 



100 - 



o 



< 50 - a.0.135 

o - -X 



-50 - 



0.121 



2 



5 



10 



20 30 50 



100 



200 



500 



nig- [GeV] 

FIG. 7: Values of ^-Xtot vs m g are shown for for different values of a s (Mz) for the 2004 study (blue 
dashed) and our new one (red solid). 

a variable a s {Mz)- In this case, the constraints from both the hadronic scattering and 
direct measurements of a s [Mz) are relevant. The most meaningful log-likelihood term is 
Xtot = Xh.s. + Xa s - If Oi s {M z ) varies, the hadronic scattering data on their own, including the 
HERA-1 and Tevatron jet data sets, are compatible with practically any gluino mass. The 
contribution to x 2 from the hadron scattering experiments, xl s > stays approximately the 
same as in the SM case, or improves slightly, as gluinos with masses of 10, 20, and 50 GeV 
are introduced. 

This agreement with the hadronic scattering data, hardly affected by the gluino mass, 
results from compensation between modifications in the shape of the gluon PDF and an 
increase in the preferred a s (Mz), which grows from 0.118 in the SM case to 0.132 for m g = 10 
GeV. In contrast, the total likelihood function for the fit to the hadronic scattering data and 
a s values, introduced as Xtot m Eq. (llOjl . varies considerably as a function of the gluino 
mass. Our assumed high-Q constraint of a s [Mz) = 0.123 ± 0.004 is slightly higher than 
a s (M z ) = 0.118 obtained by the SM evolution from a s (Q = 5 GeV) = 0.213 ± 0.002 in 
Eq. ([2]). This enhanced value of a s [Mz) would favor a slower QCD evolution above the 
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gluino mass threshold at about 50 GeV, cf. Fig. [2j Consequently, Xtot * s smaller at rrig = 50 
GeV than in the SM case, with the difference dependent on the choice of the high-Q value of 
a s (M z ). Specifically, we observe that XsusY( m g) ~ Xctio can be as small as —50, if we take 
a s {Mz) = 0.123 ± 0.004, but this difference decreases if a smaller a s (Mz) is used for the 
high-Q constraint. For lower gluino masses of 10 or 20 GeV, a s (Mz) increases and eventually 
is incompatible with the direct constraints. 



A. Ax 2 as a function of gluino mass 

The behavior of Ax 2 in the whole range of gluino masses is illustrated by Fig. [6] for a 
fixed a s [Mz) — 0.118, and by Fig. [7] for a floating a s (M z )- The quantitative likelihood of 
a given mass m g is specified by Ax 2 = X 2 { m g) ~ Xctic the difference from the x 2 value 
obtained in the CT10 SM fit. Values of Ax 2 in excess of 100 units disfavor an assumed m g 

391 ] . while a negative Ax 2 indicates a preference for this m g . 
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at about 90% C.L. Ref. 
Variations in Ax 2 with a magnitude below 100 units can result from a variety of sources and 
are generally viewed as not significant enough to warrant strong conclusions. 

In Fig. [61 two curves are shown for Ax\ s , the difference between the log-likelihoods in 
the fits performed in the SM+SUSY and SM scenarios for a s [Mz) = 0.118. Here Ax\ s 
is computed from the hadronic scattering contribution only, x\ s ■ The blue (dashed) curve 



represents the 2004 analysis 23]. The red (solid) curve is obtained in the present study, 
resulting in a tighter lower bound on m g . The left branch of the 2010 curve intercepts 
the Ax\ s = 100 line at m g « 25 GeV. The 2004 curve allows for 15 GeV gluinos and 
has a broader valley with respect to the 2010 one. This figure shows the improvements in 
the constraints from the present study, reflecting the inclusion of the latest precise data and 
technical advances in the the CTEQ analysis since the 2004 publication, including treatment 
of correlated systematic uncertainties and normalization uncertainties. 

Figure [7J illustrates the fits with a variable a s (Mz)- Two curves are shown for Ax 2 ot , 
the difference between the log-likelihoods in the fits performed in the SM+SUSY and SM 
scenarios in 2004 (blue dashed line) and 2010 (red solid line). Best-fit values of a s (Mz) for 
some gluino masses are indicated by numerical labels near each curve. In this figure, Ax 2 ot 
is computed from the total function x\ot- ^ includes the direct constraints on a s (Q) in the 
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current study and does not include the a s constraint in the 2004 fit. 

The figure emphasizes our earlier observation that the direct a s constraints improve the 
constraining power of the global analysis. At nig — > oo, the fit converges to the pure QCD 
value and a s (Mz) ~ 0.119. According to the A% 2 < 100 test, gluinos lighter than 15 GeV 
are disfavored for all a s [Mz)- Gluinos in the mass range 15 to 50 GeV are allowed if a s (Mz) 
takes a value in the range 0.121 to 0.131. Gluinos heavier than 50 GeV are allowed for 
practically any a s [Mz) value. By contrast, the 2004 curve exhibits only a shallow minimum 
around 5 to 6 GeV, and it is relatively flat as compared to the 2010 curve. The 2004 curve 
does not establish pronounced lower bounds on m g , for a free a s [Mz)- 

The 2010 curve in Fig. [7] exhibits an intriguing minimum for a gluino of about 50 GeV, 
corresponding to a s [Mz) of 0.121. Other that noting it, we choose not to base conclusions 
on this minimum for two reasons. First, from the point of view of the fit itself, given its 
initial inputs, we adhere to statement that only values of |Ax 2 | in excess of 100 units are 
considered significant. Second, the depth of this minimum is a reflection of the value of the 
input constraint a s [Mz) = 0.123 ±0.004. The dip grows deeper (becomes more shallow) if a 
larger (smaller) value of the direct constraint is taken at Mz- For example, gluinos with mass 
50 GeV would be disfavored if the direct constraint a s (M_z) < 0.118 were taken, compatible 



with some existing analyses of LEP data in pure QCD 2l|, |22j. Stronger conclusions on m- g 



await an independent reduction in the uncertainties on a s {Mz) 



B. Comparison with Tevatron jet cross sections 

Table HT1 indicates that the hadronic scattering data, including the combined HERA-1 and 
Tevatron jet cross sections, may still allow contributions from fairly light gluinos. This ob- 
servation is somewhat counterintuitive with regard to the precise Tevatron jet cross sections, 
which could be expected to be sensitive to non-SM contributions in the strong interaction 
sector. SUSY degrees of freedom introduce new subprocesses in the jet cross sections, such 
as gg — > gg and qq — > gg. The change in a s (Q) and the alteration of the gluon and quark 
PDFs also influence the jet rate. However, jet cross section measurements are affected by 
systematic effects that dominate over statistical uncertainties, notably by the uncertainty on 
jet energy scale and jet energy resolution. Correlated systematic shifts in the Tevatron jet 
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CDF Rim-2 jet production, Vs= 1.96 TeV 
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FIG. 8: Comparison of theoretical predictions for single-inclusive jet cross sections with experimen- 
tal data from CDF Run-2 for two bins in jet rapidity y. 
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D0 Rim-2 inclusive jet production, Vs= 1.96 TeV 
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FIG. 9: Same as Fig. [8l for two bins in jet rapidity y in the DO Run-2 jet cross section measurement. 



21 



data must be taken into account when comparisons are made to theory predictions [44|]. In 
our study, systematic uncertainties in the jet data limit the strength of our conclusions. 

These observations are illustrated by plots of the CDF Run-2 and DO Run-2 data vs. 
theory in Figs. Eland ED Our results are computed with a floating a s . As reference values, 
we use SM cross sections computed with the CT10 PDFs. Differences from the SM cross 
section are presented as 

((Jj - <T C Tlo)/o"CT10, (15) 

where <jj are the SM+SUSY cross sections computed for gluino masses of 10, 20, and 50 
GeV. The values of the jet transverse momentum px are displayed along the horizontal axis. 
Two bins in the rapidity variable y are shown for each experiment; the behavior in the rest 
of the bins is similar. 

The lower (red) error bars represent the unshifted data. The upper (blue) error bars show 
the data that are shifted by their systematic uncertainty so as to maximize the agreement 
with theory for m g = 10 GeV. Without the correlated shifts, the data would disfavor the 
light gluinos with a mass of 10 GeV. The perspective changes significantly if systematic shifts 
are allowed: the line representing m g = 10 GeV now lies completely inside the error bars. 
Similarly, if m g is equal to 20 or 50 GeV, the effective shifts of the data change to achieve 
acceptable agreement with the theory curve for this mass. 4 

The systematic uncertainties make it difficult to disfavor the light gluinos solely on the 
basis of the Tevatron Run-2 jet data. The figures show that the gluino contributions affect 
the whole pr range, as a result of the momentum sum rule and other connections between the 
PDFs of different flavors and at different (x, Q) values. Modifications in the jet cross sections 
due to "new physics" associated with the gluinos cannot be isolated to a specific px interval, 



contrary to the assumptions made in some experimental studies of jet cross sections 30] . 



VI. CROSS SECTIONS AT THE LHC 



The possible existence of color-octet fermions with masses in the range 30 to 100 GeV, 
allowed by hadronic data according to our analysis, raises the prospects for their detection in 



4 



The extent of plausible systematic shifts is determined by matrices of correlated systematic errors, provided 



by both Tevatron collaborations and implemented in the CT09 44[ and CT10 analyses 25 [ 
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the extended range of transverse momenta at the LHC As explained in early sections of this 
paper, these new fermions modify QCD parameters, primarily the QCD coupling a s (Mz) 
and the gluon and sea-quark PDFs. Precise studies of cross sections at LHC energies thus 
have the potential to reveal differences from pure SM QCD, such as the presence of color- 
octet fermions, provided the LHC measurements are supplemented by a robust program to 
reduce uncertainties in a s , PDFs, and other SM parameters, which may otherwise reduce 
sensitivity of the LHC observables to the gluino contributions. 

Compare, for example, single-inclusive jet cross sections at the LHC energies a/s=7 and 
s/s=14 TeV, computed at NLO with the EKS code 45|, |46j in the pure SM case and in the 
presence of light gluinos. The CT10 asymmetric PDF error bands on the cross sections, 
normalized to the predictions based on the central CT10.00 PDF set, are also shown in 
Figs. [TOirrSI as a function of the jet's transverse momentum px, in several bins of the jet 
rapidity y. Ratios of the expectations based on the SM+SUSY PDFs for m g = 20 and 50 
GeV to their counterparts based on the CT10.00 set are shown as the dashed and dot-dashed 
lines, respectively. 

In Figs. ITOlandim these ratios are computed with a s (Mz) = 0.118 assumed in all PDFs 
and cross sections. In this case, the SM+SUSY curves lie outside the respective CT10 PDF 
uncertainty bands for some pt, suggesting that the SM and SM+SUSY scenarios can be 
distinguished, if sufficient experimental accuracy is achieved. On the other hand, if a s (Mz) 
takes the values of 0.126 and 0.121 that are preferred in the SM+SUSY fits with m g = 20 
and 50 GeV, respectively, then the SM+SUSY curves lie within the CT10 PDF error bands, 
as shown in Figs. fT2l and [T3l In this case, discrimination of the SM and the SM+SUSY cases 
is more challenging, as reduction of the experimental uncertainty below the current PDF 
uncertainty would be necessary. 

For the inclusive jet cross sections to provide a good discrimination between the SM and 
SM+SUSY scenarios, the uncertainties on both a s and PDFs must be reduced below the 
current values. NNLO contributions to SM processes and NLO gluino contributions must 
also be implemented in both the PDFs and jet cross sections. 

A different approach to detecting the presence of new colored states could be based 
on the expectation that QCD radiation off a heavy colored object differs from that from 
massless partons that dominate the inclusive cross sections. It may be possible to identify 
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jets containing gluinos by studying distributions in the jet mass or other jet shapes. The 
distribution in the jet mass produced by conventional QCD radiation decreases smoothly 
as the jet mass increases. Decays of gluinos would result in jets whose mass distributions 
peak at trig, and gluino jet contributions could be identifiable above the continuous SM 
background in the distributions in the jet mass or related observables, using methods being 



developed 



47 



50]. 
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FIG. 10: Ratios of single-inclusive jet cross sections at y/s = 7 TeV, obtained from the central PDF 
set of CT10, CT10.00, and the SM+SUSY PDFs for gluino masses m~ g = 20 GeV (dashed line) and 
50 GeV (dot-dashed line). The asymmetric PDF uncertainty of the CT10 set is also shown as a 
filled band. The SM+SUSY PDFs are obtained under the assumption of a s {Mz) = 0.118 for both 
sets. 
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FIG. 11: Same as Fig. [TDJ for y'i = 14 TeV. 
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LHC 7 TeV. Fit with a floating oc s (M z ) 
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FIG. 12: Same as Fig. QUI but with a s {M z ) = 0.126 in the SM+SUSY calculation with m- g 
20 GeV, and a s (M z ) = 0.121 in the SM+SUSY calculation with m~ g = 50 GeV. 
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LHC 14 TeV. Fit with a floating a s (M z ) 
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FIG. 13: Same as Fig. [El for = 14 TeV. 
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VII. SUMMARY AND DISCUSSION 



In this paper, we explore modifications in QCD scattering cross sections introduced by 
color-octet Majorana fermions in supersymmetry (gluinos) and other popular extensions of 
the standard model. Their influence must be included in the evolution of the strong coupling 
strength and the parton distribution functions, especially if these fermions have mass below 
100 GeV (possible in the absence of model-specific assumptions). In addition to modifying 
the evolution of ot s (Q) and the PDFs of the SM quarks and gluons, a relatively light gluino 
also introduces new production channels such as gg — > gg in the inclusive jet production 
case. In this context, hadronic scattering data included in global PDF analyses can provide 
model-independent constraints on the color-octet particles. 

We examine the values of x 2 obtained from our global fits as a function of the gluino 
mass frig. By analyzing a combination of the latest HERA and Tevatron data on hadronic 
scattering, and world measurements of the QCD coupling at Q < 10 GeV and Q = Mz, 
we conclude that gluinos must be heavier than 25 GeV at 90% C.L., if a s (M z ) = 0.118, 
and heavier than 15 GeV if a s (Mz) is arbitrary. These constraints supersede the 2004 
study based on the CTEQ6 data set, in which we found a lower limit on the gluino mass 
of rrig > 12 GeV for a s (Mz) = 0.118, and no limit if a s (Mz) is arbitrary 23[. These new 
bounds are comparable to the gluino mass bounds m g > 26.9 and 51 GeV obtained from 



the analysis of event shapes in e + e~ hadroproduction at LEP [l6j, [17J. Our constraints on 
m g are obtained from the analysis of inclusive QCD observables and are not affected by 
theoretical uncertainties of the kind that arise in the determination of a s [Mz) from the LEP 



data 



18H22j and LEP event shapes. 



The changes in a s (Mz) and in the PDFs of standard model partons must be taken into 
consideration when QCD tests are made with LHC data. The high energy of the LHC and 
the extended range in jet transverse momentum offers hope that BSM deviations from pure 
QCD will show up in inclusive jet cross sections. As discussed in our comparisons with 
Tevatron jet data, it will be critical to control experimental uncertainties on the jet energy 
scale and jet energy resolution. Gluino contributions and adjustments in the SM parameters 
tend to offset one another. The power of precise measurements of the LHC single-inclusive 
jet cross sections will be enhanced provided that a s (Mz) and the PDFs for gluons and quarks 
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are constrained more tightly than now by measurements in other channels. 

For the purpose of studying jet properties in detail, we provide routines to interface 
with the SM+gluino PDFs. These are linked from the CTEQ webpage at cteq.org. We 



also note that the MadGraph/MadEvent programs 51] provide a mechanism to incorporate 
SUSY PDFs in the initial state; information for using this interface is also provided on the 
webpage. 
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Appendix A: Modification of the Strong Coupling 

The running of a s (Q) must be matched to the individual PDF set with the appropriate 
mass thresholds. The expansion of the evolution equation for a s (Q) 



n=0 



a 2 a 3 
2ir 2 ci 7r / 



(Al) 
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can be solved perturbatively. It takes the form 54] 



a s (Q) 



An 



/So In 



ft ln[ln(Q 2 /A 2 )] 
ft 2 ln(Q 2 /A 2 ) 



+ 



ft 2 



(A2) 



ft 4 ln 2 (Q 2 /A 2 ) 

The beta functions, ft and ft depend on the number of active fermions and bosons. When 



supersymmetric particles are included 

ft = 11 
and 

38 



55], the first two coefficients in Eq. flA2j) are 

2 n 1 

~n f - 2n g - -rif, 



11 



13 



ft = 102 - —n f - A8n g - —n f ~ + —n g n f , 



(A3) 



where nf is the number of quark flavors, n g is the number of gluinos, and nj is the number 
of squark flavors. As the evolution proceeds across mass thresholds, these numbers and, 
consequently a s , must be adjusted. 



Appendix B: Gluino contributions to the single-inclusive jet cross section 

The leading-order cross section for inclusive (di)jet production, HiH 2 —> j{pz)j{p4)X, 
expressed in terms of the transverse momentum p T and rapidities y 3 , y 4 of the jets, is 

-j — -j — -j — = 4 PT y^ x i x 2 fJ-jr) /g 2 -» 3 (x 2 , Hp) ^ \M PlP2 ^ P3P4 \ 2 , (Bl) 

i/3 i/4 spin 

where x\ = mx j 1 \fs{e m + e m ), and x 2 = ttlt / ' \fs{e~ m + e~ VA ) are the parton momentum 
fractions; m T = p T + m^ is the gluino's transverse mass; and i/s is the collider center-of-mass 
energy. In our analysis, scattering amplitudes for subprocesses with gluino pair production, 
99 ~^ 99 an d qq — > 99-, are included with full dependence on gluino mass m g . Scattering 
amplitudes for the other LO subprocesses (with at least one initial-state gluino) are evaluated 



in the m g = approximation, in accord with the S-ACOT factorization scheme 34j, 135 ]. 



SUSY contributions with full mass dependence can be found in the literature (e.g., in 



56| and 37|]), but they are presented here in a consistent notation for completeness. In 
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terms of the usual parton-level Mandelstam variables, s = (pi + P2) 2 , t = (p± — P3) 2 , and 
u = (pi — P4) 2 , the square of the amplitude for qq — > gg is 



\-M- q q^gg\ 2 — g 



sm\ 4(m~ - i) 2 3 (sm~ + (rnj - t) 2 ) 4(m 2 g - uf 

3(m 2 - - t)(m? - u) + 3(m 2 - - t) 2 s(m| - t) + 3(m? - u) 2 

3 (sm? + (m| - w) 2 ) 3 (2sm? + (m? - t) 2 + (m| - w) 2 ) 



s(m| — u) s 2 



(B2) 



Here rriq is the mass of the squark, and the pre-factor 8/9 is a color factor. We report the 
expression with all the fermion mass dependence, but in our computations we have taken 



the limit rriq — > 00. 



The square of the amplitude for gg — > gg is 

„ 9m~ 9m? 21um\ 

\ JVi gg^gg\ ~ -a »*or* ?.\ ^ 



1 4s 2 (t - mf ) 4s 2 (u - mf 4s 2 (t - mf 
45mi 2limi 45m^ 

j 



2s (i — mf 4s 2 (u — mf 2s (« — mf) 

21 m A h 9mi 81m^ 81m£ 

j 



(i-mf(u-mf s 2 {i-ml) 2 (u - mf)' 
27u 2 m 2 g 9tm~ Adiim 2 , 9um~ 9mf 



4s 2 (i - mf) s 2 2s (i - mf) s 



2 



s 



27t 2 m? 45tm? tf? 9u 2 9tu 

+ TTTT + TT7T + 



4s 2 (u-m 2 -) 2s(u-mf U 2 As 2 2s 2 
9u 3 9*3 



4s 2 (£ - mf ) 4s 2 (m - mf) 
Appendix C: Parton Luminosities 



(B3) 



Parton-parton luminosity functions portray the relative size of various partonic contribu- 
tions. The parton luminosity is defined as a convolution integral of the PDFs fi(£,Q) for 
two incoming partons = g, g,u,d, s, ...): 



dr 

where r = s/s. Here s is the square of the center of mass energy in the incident parton-parton 

system. In terms of this luminosity, the production cross section for a specific reaction is 

dCij(r,Q) 
dr 



£ / dra^r) dL ^K (CI) 
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m~=50 GeV, Q= 100 GeV 
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FIG. 14: Parton-parton luminosity r dCij(r,Q) / 'dr vs. \fr for m g = 50 GeV at Q =100 and 300 
GeV. 

The sum is over the initial-state parton flavors i and j, and CTy(r) is the partonic cross 
section for the subprocess initiated by partons 

The luminosities for some flavor combinations are shown in Fig. [14] for m g = 50 GeV. At 
Q = 100 GeV all gluino luminosities are smaller than the SM luminosities, but they grow 
in magnitude as Q increases. The gluon-gluino luminosity is roughly the same as the gluon- 
bottom quark luminosity, as would be expected from the momentum fractions presented in 
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Table I. At Q = 300 GeV the g®g contribution is comparable to that of the ordinary quarks. 
The g <S> g combination is smaller than s <8> g throughout the x range for Q = 100 GeV. At 
Q = 300 GeV, the evolution of the gluino is enhanced, and g®g exceeds various SM pairings 
for x > 0.1. 
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